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ABSTRACT 


Dynamic instability of planing craft on calm water, porpoising, 
is a phenomenon which has not been properly understood. Empirical 
relations are available for predicting the regime of stability. The re- 
lations, when compared, lead to conflicting design requirements to in- 
crease stability. 


It is therefore desirable to develope a theoretical approach to 
the problem so that the effects of beam, deadrise angle, etc. on 
stability can be studied. 


The results of the investigation imply that a decrease in dead- 
rise angle, a decrease in beam and an increase in distance from LCG 
to transom result in an increase in stability. Changes in shaft angle 
and vertical height of the center of gravity and moment of inertia have 
very little effect on the stability of a boat while it is planing. However 
further investigation is required to verify these results. 


In conjunction with this paper, a computer program was written 
which can be used in the design of planing craft to predict boat attitude, 
wetted surface area, drag and effective horse power. This program 
will be available for use in the XIII Department library. 


Thesis Supervisor: Philip Mandel 


Title: Professor of Naval Architecture 
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NOMENCLATURE 
Standard 
Symbol Definition 
b Beam 
CH Lift coefficient for prismatic surface 
a Lift coefficient for zero deadrise surface 
D, Drag force tis” 
F Froude number °// eV 
KG Height of center of gravity above base 
line (ft) 
} Non-dimensionalizing length (ft) 
Location of center of pressure forward 
cP of transom (ft) 
lag Location of center of gravity forward of 
transom (ft) 
lrn Mean wetted length (ft) 
L, Length of wetted chine (ft) 
Ly. Length of wetted keel (ft) 
mr* Boat pitching moment of inertia about CG 
mX, Added inertia effect about Y-axis 
N Normal force (lbs) 
T Thrust 
Mean velocity of flow past bottom 
Angle of keel above horizontal (deg) 
Non-dimensional velocity (fps) 
Wetted surface (ft*) 
Average deadrise angle (degrees) 
A ‘Displacement (Ibs) 
Shaft angle (degrees) 
Pitch angle (degrees) 
Ratio of mean wetted length to beam 


>» 


Note: Reference (10) uses this definition 
while reference (11) uses its reciprocal. 
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CLB 
CLO 
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Standard Program 


Symbol _ Definition Symbol 
p Mass density RHO 

7 Trira angle TAU 
(m - Zz) Vertical force per unit vertical acceleration Al 

zy Vertical force per unit vertical velocity Bl 

Z Vertical force per unit vertical displacement Cl 


S, + mx.) Vertical force per unit angular acceleration D1 


a, + UoZ.) Vertical force per unit angular velocity El 
(Z, + Uy Zz) Vertical force per unit angular displacement Gl 
a, - m,) Pitching moment per unit angular acceleration A2 
ms, + Us, M.) Pitching moment per unit angular velocity B2 


(My, 2 Oy M.) Pitching moment per unit angular displacement C2 


(Me + mX,) Pitching moment per unit vertical acceleration D2 


MY Pitching moment per unit vertical velocity E2 

My, Pitching moment per unit vertical 
displacement G2 
A1/(0.5. RHO. °) All 
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I. INTRODUCTION 


The unstable motions of planing craft have been under study for 
many years and have been the subject of much literature (see biblio- 
graphy). The ability to be able to predict the stability characteristics 
of a particular hull in the early stages of design is of importance to 
naval architects. A knowledge of the effects of variables such as beam, 
deadrise angle, etc. on stability would permit intelligent corrective 


action to be taken to increase the dynamic stability of existing craft. 


The problem of planing craft stability involves many variables 
and empi rical relations between some of the design variables have 


been developed to predict dynamic stability. 


Two formulae recently developed emperically from expIirimental 
data, (2) and (12), result in conflicting design requirements to increase 
stability (see Appendix C). It is therefore desireable to develope a 
theoretical approach to the problem so that the effects of design var- 


iables can be determined independently of experimental data. 


Perring (10) attempted a theoretical approach. His lack of 
success can be attributed to a number of causes. The foremost of 
these being lack of sufficient experimental and theoretical information 
to predict the stability derivatives accurately and the ommission of 


important terms. 
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YI. THEORY 


A planing hull, as a rigid body, has six degrees of freedom. This 
study treats the boat as a two degree of freedom system by investigat- 
ing what is considered to be the most important motions, heave in the 


Z-direction and pitch about the Y-axis. 


The equations of motion are nonlinear. To facilitate the solu- 


tion of these equations it is necessary to linearize them. 


Linearized equations of motion for ships have been developed 
by Abkowitz (1), Korvin-Korvosky (7) and others. Those of Abkowitz 
are most complete. If the coefficients, i.e. stability derivatives, are 
substituted into these equations and then the result transformed into 
the frequency domain, it should be possible tc evaluate the stability 
by using the Routh criterion, (5). 


The method used to predict stability or lack thereof proceeds 
as follows: 
1. The stability derivatives are determined, see Appendix A. 
2. The stability derivatives are substituted into the linearized equations 
for ship motion. (1) 
3, The resulting equations are transformed by substitutions of the form 
st st 
z=Zmaxe and @= Omax e 


4. An equation in S is obtained. 


The Routh discriminant is evaluated for the fourth order equation 


in S, see Appendix A. 
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II. DESCRIPTION OF WORK ACCOMPLISHED 


The firststep toward a solution was to determine the stability 
derivatives and combine them to form the coefficients of the linearized 
equations of motion, see Appendix A. The coefficients were then non- 


dimensionalized using beam as the non-dimensionalizing length (12). 


A computer program was written to solve for the Routh discrim- 
inants, see Appendix H, using as input data the results from a series 
of tests run at the David Taylor Model Basin (2), see Table 3. The 
resulting discriminants were then plotted against speed showing a 
consistency in the directions of the paths, however there was no ob- 


vious difference between stable and unstable! boats, see Figure 9. 


At this point, an attempt was made to determine the roots of the 
fourth order stability equation to examine their loci using a computer 
program from the MIT "SHARE" library. (SHARE No. 1514 RTSCH). 
RTSCH proved to be unsatisfactory. The answers obtained from this 
program are seriously in error even for the simplest of input equa- 


tions. 


It was then decided to vary in turn what seemed to be the most 
important variables: DIFB1, DIFC1, B2, Dl, D2, E11, Gi, and G2. 
This was done to determine the effect of changes in their magnitude 
on the Routh discriminant. From Figure 10 it can be seen that vary- 
ing G2 roughly grouped the stable and unstable boats with the unstable 
group centered about G2 equal to 0.38 x G2 at the point of zero Routh 
discriminant. The program was then run with G2 equal to 0.38 x G2 
go that the loci of the discriminants could be examined. The results 


are shown in Figure 11. Based on these results it was decided to in- 


1, Unstable boats, as refered to in this paper, are those which porpoised 
ata Fy less than 6.0. Stable boats are those which had not porpoised 
' before maximum test speed (2) was attained (Fv = 6.0). 
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vestigate G2 further. Since the term G2 is made up of two parts, force 
Aerm) + arm x Morse) it was decided to vary these two parts _ 
independently to see if better correlation could be achieved in either 
of the two groups. 7 Correlation was not improved, see Figures 12 
through 17. However it was observed that the discriminants for the 
stable boats changed very little with changes in G2, or its parts. On 
the other hand the discriminants of the unstable boats changed a great 
deal with changes in G2. This lead to the conclusion that there must be 
a term in the coefficients of the characteristic equation which over- 
powered G2 when the boat was stable but was of the same magnitude 


as G2 when the boat was unstable. 


Based on information obtained thus far, it was decided to investigate 
each term of the Routh discriminant to see which of the coefficients were 
controlling for the stable and unstable boats. Values of the coefficients 
obtained from program 1 in Appendix H were inserted into each term 
manually and inspection of the results was unfruitful. No obvious dif- 
ference could be detected between stable and unstable boats. It was 


concluded that the interaction was much more subtle. 


A closer inspection of each term indicated that the coefficients e, 
and M. , Dl and D2, may interact with important effects and this became 
the final step in the investigation of the coefficients. The results are 
shown in Figures 18 through 20. At this point the investigation of the 


coefficients of the equations of motion was terminated because of time 





1. At this point it was necessary to reduce the number of plots to 
three stable and three unstable boats in order both to simplify the plots 
and to make more efficient use of computer time. The unstable boats 
were selected by choosing two which had axis intercepts fairly close 
together and third whose intercept was remote from these (models 
4665-3, 4666-13, and 4668-9 in Figure 10). 
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limitations. 


Concurrently with the above work, a program was written to solve 
for all the hydrodynamic performance characteristics of the planing 
hull: planing angle, wetted surface,resistance, power requirements, 
and stability. This program uses design parameters as inputs and 
provides an easily readable output. Development of the criterion for 
the equilibrium planing condition is shown in Appendix E and details 
of the program are contained in Appendix H. In order to facilitate the 
writing of this program, it was necessary to determine an expression 
for mean bottom velocity based on an emperical function of deadrise 


angle, Appendix D. 


The coefficients of the equations of motion and the’ Routh discriminant 
based on the computed planing conditions were compared with those 
based on experimental data. The result of this comparison is shown 
in Table 2. 


The program was then run, for model 4668-9 with G2 equal to 0.38 
x G2, varying BETAI, EPSILI, VCG, BEAM, CG, and YI inturn. The 
results were then plotted, Figure 21,80 that a comparison of the re- 


lative effects of the variables could be made. 
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IV. DISCUSSION OF RESULTS 


The plots of the discriminant versus speed obtained from the first 
calculations, Figure 9, are disappointing. According to these results 
most models were stable throughout the entire range of speeds in- 
vestigated, contrary to the experimental results. The lack of agree- 


ment could be caused by one, or both, of the following: 


(a) incorrect formulation of one or more of the stability derivatives 


(b) neglecting the cross coupling effects of longitudinal motion. 


Perring (10) indicated that inclusion of the longitudinal motion cross 
coupling effects had negligible effect on the outcome of his solution to 
the problem. It is possible that in the present, more refined solution, 


the magnitude of this cross coupling may become relevant. 


The investigation of the effects of varying the magnitudes of the 
stability derivatives indicates that a solution to the problem may lie 
in this area. Although variations in Ze i, and (Z + Up Zs DIFB1, 
DIFC1i, and Gl, failed to yield any evidence of consistent influence, 
Figure 10 shows that variation of G2 produced a fairly consistent dif- 
ference between stable and unstable boats. It is true that there is con- 
siderable scatter of the axis intercepts within each group, but there is 
an undeniable consistency in the grouping. It is also of interest to note 
there is much less scatter in the stable group than in the unstable group. 
The grouping of the stable boats cannot be attributed to their equal Froude 
number. An examination of the unstable boat grouping, Figure 10, in- 
dicates that models 4665-3 and 4668-9 intercept the axis at the same 
point with Fy of 3.24 and 5.03 respectively whereas model 4666-17, 
which has a Fy. of 5.01 (essentially equal to that of model 4668-9), has 


an intersection remote from the preceeding two. 


The investigation of the loci of the discriminants with G2 equal to 
0.38 x G2, Figure 11, show that the original curves, shown in Figure 9, 
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now bend down toward negative values of discriminants as speed in- 
creases. Some of the unstable models have loci consistent with ex- 
perimental data, i.e. the Routh discriminant heads toward negative 
values as the porpoising speed is approached however, the boats which 
were stable throughout their test range have ioci lying completely below 
the axis (negative discriminants indicate instability). This indicates 
that simply multiplying G2 by a single factor does not produce reliable 


Routh discriminants. 


The investigation of the effects of varying the two components of G2 
separately, Figures 12 through 17, although not producing better corre- 
lation within the stable or unstable groups does point out the small ef- 
fect that these variations have on the intercepts of the stable group 
compared to their effect on the unstable boats. This information, as 
it stands, indicates that other terms in the equations of motion are 
more powerful at stable speeds but that, at the porpoising speed, G2 


is a powerful term. 


The results of the investigation of the simultaneous variations of Dl 
and D2, see Figures 18 through 20, point the way to what may be a 
valuable area for further study. The first useful bit of information ob- 
tained is the fact that oY D1, has very little effect on the magnitude 
of the Routh discriminant and that M., D2, has a large effect. The 
most important result of this investigation is the fact that the axis in- 
tercepts of models 4666-13 and 4668-9 have been reversed in their 
relative positions from what they were when the coefficient G2 was 
varied. This means that a simultaneous variation of D2 and G2 may 
cause these two extreme boats to cross the axis at the same point and 


thereby correlate the unstable group. 


Correlating the results in this manner does not really solve the 
problem. The stability indicator, Routh discriminant, is not reliable, 


as the discussion of Figure 11, Appendix F has shown. Further work is 
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required to produce better agreement between the intercepts of the 
discriminants and experimental data. An experimental investigation 

of the individual stability derivatives for comparison with the theoreti- 
cally developed derivatives would be helpful in locating the terms of the 


equations of motion which need to be reevaluated. 


The program, which solves for the hydrodynamic performance 
characteristics of planing boats, yields information of importance to 


design. 


Starting with attitude, wetted keel length, wetted chine length, and 
drag, it can be seen, Table 1, that there is good agreement between 
theory and experiment for boat attitude and drag. The theoretical values 
of wetted keel length and wetted chine length are larger than the ex- 


perimental values. 


A plot of WKEEL - WCHINE vs TRIM comparing theory (12) and 
values calculated from the experimental results of (2), Figure 22, in- 
dicates that the mean line of data points lies above the theoretical line 
for this group of boats, Table 3. The largest errors occur at the largest 


angles of attack. 


The expression developed for mean bottom velocity, Appendix D, 
yields results which compare very accurately with graphs shown in 


Figure 7. 


The comparison of derivatives calculated directly from the program 
shows good agreement with the exception of C2 and E2, see Table 2. 
This error was most likely caused by the difference in actual and 


calculated wetted length. 


The results of the variation of BETAI, EPSILI, BEAM, YI and CG, 
shown on Figure 21, can not be conclusive because of inconsistencies 
. which have been found in the discriminant. However Figure 21 does show 


that the beam, deadrise angle and longitudinal position of the center of 
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gravity have the largest effect on stability. The implication of Figure 21 
that a decrease in beam increases stability agrees with (12). The in- 
ference that an increase in deadrise angle results in a decrease in 
stability is, at first glance, distressing. However an inspection of 
Figure 16 of reference (12) indicates that this may well be the case 

for the boat examined. An increase in deadrise angle results in an in- 
crease in trim and a decrease in the lift coefficient. Both of the latter 
effects are destabilizing. If the destabilizing influence caused by the 
change of trim and of the lift coefficient is greater than the stabilizing 
effect of the increase in deadrise angle, the boat will be destabilized. 


The results indicate that moving the center of gravity forward in- 
creases the range of stability. This forward movement results ina 
decrease in trim and it is known that a decrease in trim results in an 


increase in the range of stability. 
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V. CONCLUSIONS 


Unfortunately this thesis does not go far enough to solve the problem 
of predicting the stability characteristics of a planing craft in the de- 
sign stage. The relative magnitudes of the force and moment derivatives 
which make up the coefficients of the equations of motion give rise to 


inconsistencies in predicting stability which have not been resolved. 


This study does indicate however that variations in the following 
quantities have a minor effect on the magnitude of the stability in- 


dicator (Routh discriminant): 


(a) change in lift coefficient x area with respect to vertical velocity 
(DIFB1) 

(b) change in lift coefficient x area with respect to vertical position 
(DIFC 1) 

(c) change in vertical force with respect to the angular acceleration 
about the pitch axis ia, or D1) 

(d) change in vertical foree with respect to the angular velocity about 
the pitch axis a" + Uo. VERM or E11) 

(e) change in vertical force with respect to trim (2, + Uy Zz. or Gl) 

(f) change in pitching moment with respect to velocity in the heave 
direction aa + Uy M. or B2) 

and that the following have a major effect: 


(a) change in pitching moment with respect to vertical acceleration 
(M. or D2) 

(b) change in pitching moment with respect to vertical position 
(M, or G2). 


The computer program, developed as part of this thesis, which solves 
for the other hydrodynamic performance characteristics of planing craft 
is able to reproduce experimental results with minor limitations. The 


| expressions used in computing wetted length of chine and keel do not 
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yield accurate results. This causes an error in wetted length and 
aspect ratio which is further reflected in some of the coefficients of 


the motion equations. 


The expression developed for computing the mean bottom velocity 


yields consistently good results. 


The results of the study of the variations of design parameters 
generated by means of computer Program 2, Figure 21, show that the 
range of stability may be increased by decreasing deadrise angle, de- 
creasing beam and moving the center of gravity forward. Changes in 
the vertical center of gravity, moment of inertia and shaft angle have 


minor effects on the range of stability. 


In spite of the inconclusive results of the stability investigation, there 
is substantial indication that the stability problem can be solved. The 
solution of this problem will utilize results like those shown in Figure 
21 in conjunction with Program 2, This should provide an extremely 
useful design tool for optimizing planing craft design. For example; 
assume that it is desired to design 2 planing hull to operate at 40 knots. 


The design procedure would proceed along the following lines: 


(a)Run computer Program 2 for a number of combinations of beam, 
deadrise angle, and the longitudinal position of the center of gravity 
obtaining drag information for all combinations which yield designs 
stable to 40 knots. 

(b) From the data thus obtained, develop a family of curves for each 
deadrise angle by plotting drag versus beam for several locations 
of center of gravity. 

(c) Choose the design for minimum drag for a 40 knot planing hull from 


the curves. 
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VI._ RECOMMENDATIONS 


i. Repeat the study described herein using three degrees of freedom: 


pitch, heave, and surge. 


2. Develop a more spaare* method of predicting wetted length of 
of CAINE 
keel and wetted length’for a planing surface with deadrise angle based 


on experimental data. 


3. Make an experimental investigation of the force and moment 


derivatives for comparison with those developed theoretically. 


4. An examination of the loci of the roots of the characteristic equa- 
tion in the S-plane would produce valuable results once the inconsistencies 
in predicting stability are ironed out. An investigation of this sort would 
show the effect that variations in design parameters have on how the 

roots approach and cross the imaginary axis. This requires a more 
accurate root extraction computer program than was available to the 


authors. 
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APPENDIX A 


Equations representing force and moment equilibrium for a ship can 
be expressed in the form: (1) 
AlZ + BIZ + C1Z + DIO + BI@ + Gia + Fe __.. forces 
A28 + B20 + C260 + D2Z + E22 + G2zZ = ue’ .... moments 
For the case under consideration here, smooth water, there are no 


external excitations. Therefore re™ and ure are both equal to zero. 


By substituting the transforms ze™ and ge%t into the force and mo- 


ment equations we can obtain a characteristic equation of the form: 


AA. S*+ BB.S° + cc.3° +DD.S+ EE =0 
Where: 
AA = 1. 


_ A22,.B1l + Ail. B22 - D22, Ell - E22. D11 


—_ Al1.A22 - Dil. D22 

oc = £22.C11 + B22, Bll + All. C22 = D22. Gil - B22. Ell - G22. dil 
‘ Ail A22-Dil p22 

pp = b22. Cll + Bll. Cee - B22. Gil - G22. Hil 


A11,A22 - D11.D22 


C22.C1l - G22,Gil 


EE * 411. Aaa - D11.D22 


If the Routh criterion is applied to the characteristic equation it should 
be possible to evaluate the stability of the boat (5). The criterion indi- 
cates that a boat will be stable and nonoscillatory in the steady state if: 


BB.CC. DD - AA. DD’ - BB’. FE > 0 


Negative Routh discriminants are not meaningful for the case under 
study. Negative roots denote instability, instability implies motion, and 
motion in this case implies nonlinearity. Since the method is based on 


linearized equations, the last meaningful Routh discriminant is zero. 
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THE DERIVATION OF COEFFICIENTS FOR FORCE AND MOMENT 
EQUATIONS: 


Al 


Al and AZ will have to have terms containing added mass and added 
inertia respectively. . This is necessary because the added mass and 
added inertia may be of the same order of magnitude as that of the 
boat itself. 

Al =W+ VERM 


Bl 
Bl = Z 

W 

AZ _ 30. aZ 
Z = owe FS 
w ow ow 08 
ee. i 
aw V 
aZ 1 g eCLA 
= 9 = @ RHO. uo. a 0 








Z.=0.5 RHO uy DIFBI 











Bil =Z 

Ww 
where: 
CLA = Cy A (lift coefficient x area) 
DIFB1 = Cc (DIFB1 is the symbol used in the computer program) 
C1 
Cl= Z, 

~e2 . im) 1 > (CLA) 

*. az ~C~<“‘C CD = ¢ BHO, az 


4, = 0.5 RHO u,” DIFC1 


CizZ 
Z 





* 
See Appendix B. 
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where: 


aCLA 
OZ 





DIFC1 = computer program symbol for 
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E€=VCG. cota 


* 
See reference (9), page 420, Fig. 62A. 
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q 
FIGURE 2 
ee +t ly. Zs 
q 0 
_942,9o9w 980 04 ._ dw Bl 
oq cq° ow’ 30 oqa 
now wish to express q as an effective velocity, w. 
VCG 
sina 


v2zr.q. cos(a + @)= VCG. q(cos 8. cota - sin @) 


d(v.a) = increment of velocity x area 


= VCG.q. (cos 8. cota -sin9@)d&. BEAM 








=e 


€= CG 
V.A= -VCG*, BEAM , q.[cos ofcot a , csc?ada - Sin afcsc%a dal 
E=CG-WETL 
where: A = wetted area s WETL. BEAM 
V.A=-WETL. BEAM. q.[cos @(0.5 WETL - CG) + VCG, sin 0 | 
Z_ = B1 [cos @ (CG - 0,5 WETL) - VCG. sin 0] 


El=Z +uo. VERM 





Gl 

Gl = Z, + Uo. Z 

Lift = =. RHO. ug. CLA = Z 

Zee» 1 2 o(CLA) 

56 5 ° RHO. Uo 5 8 

Z, * >. RHO. ug DIFB1 = uy B1 

Gl = Z, + Uo Zi, = Uo Bl+u, Bl 
=Z. Uo. Bl 


x 
A2 = YI+ VERYI 
where: 
YI = pitching moment of inertia 


VERYI = added inertia, 


a 
See Appendix B for development of VERYI, 
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FIGURE 3 
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FIGURE 4 
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Sie 
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r = (VCG? - CG?)” 
WK =H/sing , 7 =sin~ (VCG/r) 


H =r.sin(7+0)-(VCG-H ) 


6=0 
. 0(arm) 6 (force) 
M, = force Sons + arm ok 


= W Seer) [(CG - CPL) cos @ - VCG. sin @].G1 


Assuming:. CPL = C,.WK 


Siar) = - CG.sin 8 - VCG.cos 9 - C,ecos 9 oom + C,-WK.sin 9 
LS =-resinze csc”? @+ VCGecot @-csc 8 


= VCGecsc @ (cot @ - csc 8) 
VCG 





o(arm) _ p ’ 7 ; 

0 = CPL [sin @ waz cot 8 (cot @ - csc 9) ] - CGesin 80 - VCG-ecos 8 
C2= My + Uo M = Mg + Uo E2 

D2 





FIGURE 5 
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dM = = » RHO.BEAM?. dé (é.cos 8 - VCG-sin 6).w 


Note: Incremental force = d (added mass. w) 
added mass term taken from reference (11), page 420, 


Figure 62A, for an increment of length. 
CG 


M. = > RHO. BEAM? (&.cos @ - VCG.sin 0) dé 
CG = ‘Waals 
2 =~ -RHO . BEAM”. WETL [cos @ (CG - .5 WETL) - VCG.sin@ ] 
D2=M. 
Ee VCG- SING 





FIGURE 6 
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= (CG - CPL) cos 8 - VCG.sin 8 ] eau 


M = [ (CG - CPL) cos 9 - VCG.sin 9 ] Bl 
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See Figure 4. 








G2 = M, 
cM 3(force. arm) 

M_ = iin. % “iain 

Z og em 4 

= force df(arm) + arm LE nae | (force) 
OZ OZ 

a(arm) — rs a (CPL) 

roe 2 roe A 


where C, is an arbitrary coefficient. 

















|. @ 
=a 6 
where H = draft at transom. 
oe.) 4B! 
32 sin® é2 sin 8 
d(arm) _ _—- cos @ =—* 
“sn rm ai sin 3) C, C, cot 9 
CPL 
- “WE cot 8 
CPL 
M, = L(CG - CPL) cos 6 - VCG sin9 } C, - W Ge cote 
G2=M 
a 
8(arm) 
Note: Computer program determines “te by incrementing variables 
1 


in the equation CPL =0.75 - :; 
5.21(-—")? + 2. 39 


developed in reference (10), page 16. 
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APPENDIX B 
Solving for Added Mass: 


Because of a lack of information about deadrise surfaces the added 
mass is calculated as for a submerged elliptic cylinder. The result 
will be divided in half and the minor axes of the ellipse will be set to 
zero. This is assumed to be a suitable approximation of the added 


mass of a deadrise hull at the water surface. 


For flow past an elliptic cylinder (page 251, (8): 


~4.,8 
Pa - z 
where: w = complex potential 


y= re a unit circle 
A= U(bcosa +iasina) 


B = 3 w (a®- b) 


For the case in question where only the vertical oscillation is being 


considered, 
w= 0 
ns wat - Ulb cosa +ia sina) 
; ~y x 
w-A> 
aw - 
| <he . 
Now: T2 - L! ie/ - ay T = kinetic energy 
4 be) | 
(c) 
= - : ioe 2 wi (the residue AA) |} 
»% : i © 21 i U2(b®c08® a + a%sin%s } 
2 
For this case a = 90°, b=0 and T = - p2xa® 


2 


but T = > M U? 


then 2 VERM = pe ¢ a? per unit width 


and VERM « = (WETL)*® (BEAM) 
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The Added Inertia: 


The problem here is to determine the added inertia for a body 


rotating about some point other than its center. 


For rotation about a point other than the center of the body the 
complex potential is given by: 

We iw Ln Z 
In this case Z, = Z because it is on the real axis. 


Z2 acos 7+ { b sin 7; 
then 


we 3+5-to Zo, (a cos 77+ i b sin 7) 


w-=AT+ BY®+iw Z,a cos 4+ bw Zo, sin f; 
dw 





- dy) 
a A + 2BY - iw Za sin 7; tbe Zq cos 1; ay 
Se @® 

SS. im 6h 

@ '* “its 

dw 7 WLo * _-_ 

a> A + 2BY + — —=" (=~ b cos 5 a sin 1) 

WW = 45 + San! Oe ) 
hls datthetins ai oben 

iw*Z, 7 





(a cos i + ib sin 7) (= b cos 7;- a sin t) 


T = - 1 pfwawe 


1 AB 
se lf + — +oER +e Bas =e as 





+ - 20 (a cos 1 +ib sin) ( bb cos 4 - asin) dS 
+ i 

(c) 

+ f -i Aw Z, (a cos 7+ib sin i) df 
(c) 

a4 - £1E Ge Sete eoeqe ib ote wy) az | 
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=F1l+F2+F3+ F4 
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in this case A = A = 0 because there is no Z translation. 
Therefore: 
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for the plate, b = 0 
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APPENDIX C. 


From Figure 16 of (12) the stability of a planing craft can be in- 


creased by increasing 


G 
. 
“> i.e Cc, . 
But C, = 0120 »2/? ri} (page 10, (12) ) 


L 
Since increasing Cy increases stability, increasing A increases stability. 


In this case A= me . This means that a decrease in b will result in an 


increase in stability. 


From Figure 20 of (2) the stability criterion is 


Cc 
= = decreases the boat 


l 1 2 
<P. ZV b 1. 





Pp 
is stabilized. Therefore increasing b stabilizes the boat. 


A comparison of the two methods makes it difficult to decide what 


effect b has on stability. 
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APPENDIX D. 


In order to determine the resistance of a planing surface it is nec- 


essary to know the mean velocity over the bottom. 


In Figure 12 of (12) a method is provided for determining this graph- 
ically. Part of an analytic solution is provided which, if completed, 
would be useful in the computation of mean velocity in a computer 


program. 


The following is a resume of the development of the complete equation. 
From (12): 


1.1 
. 0.012 ‘ 
ee Db — i £ UB) 
A 5 cos 7 
where: 


VM = mean velocity over bottom (fps) 
V = the u,, the horizontal velocity of the origin of coordinates (fps) 
T = trim angle of planing surface (degrees) 


f(8) = an undetermined function. 


By changing 7 to radians and plotting 


1 
3 
ta- (Sy *eo8* = 46) 
1.1 
.0120 7 


it was possible to find an f(f) which yielded satisfactory results. 


The final result is: 
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VM 1 0.120 T ~- ; 
= 20-vaher REET gupteaetict SSeS a ee 
¥ r z cos T ~~? 
Z 
where: 80. - 50.8 
f(B) = cos? B 


7 = trim angle in radians 


8 = deadrise angle in radians. 
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APPENDIX E. 


EQUILIBRIUM PLANING CONDITIONS. 





FIGURE 8 


Note: This development follows reference (12) in general. However, the 
final result contains terms neglected by SAVITSKY. 


Summation of forces in vertical direction: 


(1) W = N.cos @ + T- sin(@+ €) - D,.sin @ 


Summation of pitching moments: 


(2) Nec + D,- a = “Lot. = 0 


Summation of forces along keel: 


(3) T.cos «=D, + W.sin9@ 
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Combining (1) and (3): 
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ew | | f ; 
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Combining (2) and (3): 


f cos € f 


Combining (4) and (5), we obtain the general equilibrium requirement: 
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Experimental Results from (2) 


from program 2 Results 
Model No. 4668 BVA ETS 


Table 1. Comparison of calculated and experimental 


hydrodynamicyemamacterisies, 
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. Data 


For Model 4668 Test No. 9 from (2) for the 19.36 


Knot Speed. 
















Derivative As computed from 
experimental planing 


conditions 


As computed from 
calculated planing 
conditions 


Table 2. Comparison of coefficients of equations of 


motion as computed from experimental data 


and data generated from computer program 2. 








+, Siu 


maple 3. Model Data(2). All models have deadrice = 91209" 


Model | Run | Weight (lbs)! Beam (ft) | Length (ft)| CG (ft) |. (at max. 
Number ) . test speed) 


A 
[ases [os [0.07 | acess | noe [ne | ae 
Bes [ ust [ co [ser | en) 
466g | 9 | 141.80 _ 1,190 8.00 2.95 5.03 





* hd e 
Note: Unstable boats, as refered to in this paper, are those which 
porpoised at F_ less than 6.0. Stable boats are those which 
had not porpoised before maximum test speed of reference (2) 


was attained (F ~~ 6.0). 











— 








PROGRAM i 


Program 1 can be used to solve for the stability derivatives, co- 


efficients of the equations of motion, the coefficients of the character- 


istic equation and the Routh discriminant using experimental data as 


input. 


Data cards are punched in the manner indicated by READ 1 and 
1 FCRMAT where: 


WwW 


ALFAO 


CG 


C 

RT 
WK 
wc 

S 
TRIM 


YI 


RHO 

U 
BEAM 
 BETAI 
VCG 
EPSIL 


weight of boat in (lbs) 

trim angle when boat is at rest (degrees) 
longitudinal position of center of gravity forward of 
transom (ft) 

forward speed (knts) 

towing force (lbs) 

wetted length of keel (ft) 

wetted length of chine (ft) 

wetted area (ft*) 

(planing angle - a,) (degrees) is change of trim from 
the at rest position 

moment of inertia about the Y-axis. Axis taken through 
center of gravity. (lb ft sec*) 

density of water (Ib sec? /ft) 

arbitrary non dimensionalizing velocity (ft/sec) 
beam of boat (ft) 


deadrise angle (degrees) 


height of center of gravity above the keel (ft) 


shaft angle (degrees). 
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COMPUTER PROGRAM ONE 
SQUATTONS OF MOTION OF -PEANING HUES 


RD=ROUTH DISCRIMINANT 

fee 88 +.E1C=ROUTH CRITERION FACTORS 
AllsBll»sETC=NONDIMENSTONAL FORCE AND MOMENT COEFFICIENTS 
Al»sBlsETC=FORCE AND MOMENT COEFFICIENTS 


READ INPIT DATA 


666 READ1 »WsAl FAQ2CGsCoRT sWK WC oSoTRIMSYI »RHOsU »BEAMSBETAL » 
1VCGsEPSIL 
1 FORMAT (F702 F502 sF402sF 502 F5e2sF4o2gP4e2 oF 5e2sF4ecteF5e2s 
1F5039F 3000503 9F4el 9402 2F 502) 
V=C#¥1.689 
WETL=(WK+WC)/200 
VERM=0125*RHO*301416* (WETL #2) ¥BEAM 
VERY 1=00625* 40625*RHO*3 4141 6% (WETL¥*4 ) #BEAM 
CV=V/SQRTF(3262*BEAM) 
A=WETL/BEAM 
@O=WETES(0. 75-1 .0/( 52.21% (CV/ANE 2 142.39) ) 
TAU=(ALFAO+TRIM) /5702956 
BETA=BETAI/57«2956 


EPS=EP SWS ie2956 


DIFFERENTIATION OF LIFT COEPRIECTENT X AREA WITH RESPECT 


TO TAU 








See Oise 


TAU1L=TAU=-0.001 
TAU2Z=TAU+0O,001 
A=1e/A 
GLVOUTH1ls/(2.e#WETL* (CVA Z) 4 ( (WCHEEZJASINE( 2exT AU ) 
1/BEAM+1¢/3e%*(2e*WC+WK) *SINF (BETA) /COSF(BETA))) 
3 CLB1=0e5*CLVOL1 
CLSA1=(1e¢5708* A¥TAUL*(COSF ( TAU] )*¥#¥2)*(1¢0—“SINF( BETA) ) 
1/(1¢0+A)+4e.0*%(SINF(TAUL]) ¥*¥2)*(COSF( TAU] )*#*3)*COSF( BETA) / 
2o60+CUB1)*S 
CLVOL2Z=le/(2e*WETLE(CV¥#2) ) HCCI WCHX2) KSINF(2e*TAU2) 
1L/BEAM4+1 oo /32e%*(2 o*WCH+WK) *¥SINF (BETA) /COSFI(BETA)) ) 
5 CLB2=0-e5*CLVOL2 
CLSA2Z=(1e¢5708*#A*TAU2#( COSF (TAU2) #¥2)%(1e¢0-SINF(BETA) ) 
1/7(1¢04A)+4e0*(SINF(TAUZ) *#2)*(COSF(TAU2) **3)*COSF(BETA)/ 
2320+CLB23*S 
CL=1¢5708* TAUXA*COSF (TAU) #¥2* (1 e-SINF(BETA) )/(1LetA)+4e* 
1LSINF (TAU) ¥#2*COSF (TAU) **#3*®COSF (BETA) /3e+(CLVOL1+CLVOL2)/4. 
TAO=TAU1 \ 
TAB=TAU2 
CO=1 oS TOB*TAUFA#COSF (TAO) HE 2¥(Le-SINF(SETA) )/ (1 e+A) +4e* 
morNr { TAO) Fe xCOSF( TAO) #*¥3*COSFIBET,.I/73~.+(CLVOL1)/2.6 
CB=1¢5708* TAURA¥COSF (TAB) ¥¥2¥( 1 e@SINFI BETA) )/(1LetA)t+oe*® 
TSINF (TAB) ##2*COSF( TAB) **#3*CUSF(BET4A) /3.+(CLVOL2)/2.6 
DADTAU=1000 .#W/ (RHO#¥V#*2 )*(12/CB-1-/C0) 


DIFB1=(CLSA2=CLSA1)/.002+CL*¥DADTAU 


Pie eKeNr in TON OF LIFT COEPF X “AREA WITH RESPECT FO 2 








- 61 - 


DELZ=0.001 
DELA=-DELZ*BEAM/SINF (TAU) /WETL*¥*2 
DELS=BEAM*DELZ/(SINF (TAU) *COSF (BETA) ) 
DELL=DELZ/SINF (TAU) 
CLVOL3=(( (WC-DELL) **2) *SINF (TAU) /BLAM+ (2 ¢*WC+WK—3 «6 *DELL ) 
L¥SINF(BETA)/(COSF(BETA)*30))/(2e*(WETL=DELL ) #(CV**2) ) 
7 CLB3=0-25*CLVOL3 
CLSA3=(165708* (A-DELA) *TAU*(COSF ( TAU) #*2)*(1e—-SINF (BETA) ) 
1/(1e+A-DELA) +46 *(SINF (TAU) * #2) #(COSF( TAU) *#3) *COSF (BETA) /3. 
2+CLB3)*(S—-DELS) 
CLVOL4=(((WC+DELL) **2)*SINF (TAU) /BEAM+(2¢*WC+WK+36*DELL) 
IPestNF (SETA) /(COSF(SETA)*3.))/(2e*( WETL+DELL ) ¥{ CV¥%2)) 
9 CLB4=065*CLVOL4 
CLSA4=(125708% (A+DELA) *TAU*(COSF (TAU) ¥*#2)*(1e-SINFI(BETA) ) 
L/(1LetA+DELA) +4 e*(SINF( TAU) #¥2)*(COSF( TAU) **3)*COSF (BETA) /3¢ 
2+CLB4)*(S+DELS) 
DIFC1] = (CLSA4=CLSA3)/2002 


ee, 


BIfFFERENTIATION OF MOMENT WITHERESPEGT 10.2 


Cpe =e7o° (WETE=DELL)/(5,21®(GV*GeAijesO/(WETL-DELL jee eee ess 
CPL2=675-(WETL+DELL)/(5e21*(CV*#BEAM) #*2/ (WETL+DELL) ##2+2 039) 
DCPLDZ = (CPL2-CPL1)/.-002 

C1=0265*RHO#(V¥#¥2)*#DIFCIL 
DMDZ=(CG-CPL)#C1/COSF (TAU) =-“W*¥DCPLDZ/COSF (TAU) 
WTCL=0e5*RHO*S*¥V¥¥2*CL 


Al=W/32e2+VERM 








SG Dies 


Bl=0e5*RHO#V#OIFB1 

C1=0e¢5*RHO* (V¥#*¥2)*DIF CI 

DI =VERM*VCG*¥(COSF( TAU) *#(CG=.5*WETL) /VCG=SINFI( TAU) ) 

El=B1l*¥(COSF(TAU)*(CG=WETL/2.)-VCG*¥SINF( TAU) ) 

Gl=V*B1 

A2=YI+VERYI 

B2=E1*( (CG=CPL) *COSF (TAU)=$VCG#SINF( TAU) ) 

C2=G1*( (CG=CPL)*#COSF (TAU) -VCG¥SINF (TAU) )+WH#(CPL¥(SINF (TAU) -VCG* 
1COSF (TAU)/SINF (TAU) /SWK#(COSFITAU)/SSINF (TAU) @1Le/SINF( TAU) ) )=CG* 
Z2SINF (TAU) =VCG¥COSF (TAU) ) 

D2=VERM*(COSF( TAU) *(CG=WETL/2¢)-VCG¥SINF( TAU) ) 

E2=B1*( (CG=-CPL)*COSF (TAU)-VCG#SINFC(TAU)) 

G2 = DMDZ 

El = E€1 + V*¥VERM 

G1=G1+V*B) 

B2 = B2 + V*D2 


C2=C2+V*¥E2 
a. 
QUANTITY Oe5 RHO CANCELLED OUT OF ALL FOLLOWING 


AL1=A1/(BEAM*#3 ) 
B11=B1/(U*(BEAM*#2) ) 
CLI=C1/(BEAM*(U*#2) ) 
DL1=D1/(BFAM#*4) 
ELL=E1/(U*(BEAM**3 ) ) 
GL1=G1/( (BEAM*U) ¥¥2) 
A22=A2/(BEAM¥#5 ) 


B22=B2/(l *(BEAM##4) ) 
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C22=C2/( (BEAM##3 )* (U#*2) ) 
D22=D2/(BEAM##4 ) 
E22=E2/(U*(BEAM*#3) ) 
G22=G2/( (BEAM*U) **2 ) 
AA=1l-e 
BB=(A22*B11+A11*B22-D22*E119"°E22*D11)/(A11*A22-D11*D22) 
CC=(A22*C11+B22*811+A11*C22-D22*G11-E 22*E11—-G22*o ie ae 
PeA22—D11*D22 } 
BO=(B22*C11+B11*C22-E22*G11-G22*E11)/ (Al 1 *A22S0 Bie 
mee CO226CT1—G22*G611)/7(AlLIL*AZ2=D1 k*D22) 
RD=BB*¥CC#DD~AA* (DD #*2)-(8B**2)*EE 
PRINT39 
39 FORMAT(17X93HA11917X93HB11917Xs3HCl1s17Xs3HD11917X93HE11» 
117X»93HG1) ) 
PRINT 409A11s8119C11sD11%E119611 
40 FORMAT(6F20-25) 
PRINT 41 
41 FORMAT(15X93HA22915X93HB22 915X 93HCL2915X 93HD22915X93HE22 » 
5% s3HG2Z2515Xs3H RD) 
PRINT42 9A22 9822 9C222D22£229G22»RD 
42 FORMAT(7F] 6-5) 
PRINT43 
43 FORMAT(14X92HBBs14X »2HCC 914X92HDD 914Xe92HEE 91 3X» 3HRD1911X» 
15HOIFB1l:11X»s5HDIFC1) 
PRINT449BBsCCsDDsEEsRD1] »DIFB1»sDIFCI 
44 FORMAT(7F16¢5) 


PRINT112WsALFAO 9CGoVeRT WK sWCoSeTRIMsYI »RHO9UsBEAMBETAIL » 


ih GGyePSil 








11 FORMAT(UFTe2 sF5 elo F4e2sF5 02 sF502sFoe2sFoeo2 oF 5ecehGecshoecs 






LFS e323 sFZeOsF 53 oF 4el sFoe2sF5e2//) 


GO TO 666 ——>) a 


|& | )\ | 
lean 
: 
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PROGRAM 2. 


This program solves for planing conditions: TRIM, ASPECT RATIC, 
WETTED KEEL, WETTED CHINE, WETTED AREA, DRAG, DRAFT at 
the transom, MEAN WETTED LENGTH, ESTIMATED EHP and the 
STABILITY INDICATOR (Routh discriminant: positive indicators imply 
stability, negative indicators imply instability. This section of the 


program does not yield satisfactory results as yet. ) 


A listing of all iterations involved in solving for the planing conditions 
and a listing of the coefficients of the equations of motion can be ob- 
tained as indicated by comments on the first page of the program print 


out. 


This program uses as input, Ist card: 
LIST 1, LIST 2, N BOATS 
FORMAT (3 I 3). 


LIST 1 and LIST 2 are as defined on first page of the program print 
out. 
N BOATS is the number of different boats to be run. 


2nd card: 


BETAI, EPSILI, F, VCG, BEAM, CG, RHO, YI, W 
FORMAT (4 F 5.2, 5 F 10.2) 
BETAI = BETA of PROGRAM 1 
EPSILI = EPSIL of PROGRAM 1 
F is the perpendicular distance from shaft center line to CG (ft). 


All other variables are as defined on page 58. 
3rd card: 


NUMBER, IDENT 
FORMAT (1 3, 5A4) 
NUMBER = number of speed cards which are to follow 


IDENT any identifying statement or symbol not to exceed 20 spaces. 


Aaa “IAT 


SOPDAA TIOGTRA MINT - colic poleaio SOF ge los Mavporg off ; 
 ~ MTIARG SARG ASAADATTRZY SO USTTHY 1266 GaTTaW 
om has TUS GVTAWITEX RYOUTI ONTIWW YAGIA .moceaat ost 
eer srotmnibot witive “manimiTwolh Good) AUT 6 ITT rine 
ont 1 ool eu? iiideiens Gqen easier: eviteges | 
a hens oy wast tehcahdtga> Minty ne WH neem 
Moalsibn. palmelg Ail) VO) Yabripe ml Gexlrved wmodticw!! Ms te pucltedl A 
do 9d Ako alent Iw BHOieuin ortl le sirerolttec> alt tn priteif 
ache rng rgeyy sui Yo syaq seri! eo) Wo Sioemrmoo vd totasibal wm t 


































O80 PA) tuqni as been onetyond alAT 
STAog * .& Tels i TD 
ASTON 


el emmiet vale taka texto beled ag $48 2° TERI or 


_ @ . 


AEM SO adMed Taw2eTih Yo =welomen all wi STM 


-_ 


WUT ORF O° (MASc 40° 4 Auece antes 
QOL FAS) TAD 
i MARDINT 10 ao, are 
i jaahrinay ys ume. nse 
iiiiarllkes deen units tachsiiesicee ira silt af 
0S opag ao bectt4d a8 aca esidehvew vs 


Aa 


rise aoa 


(bAC 61) a 


allot wt ote Uelaw dian heegs Yo Ledeen > RSUMTM 
Lemaye OF Rpncee oF ton Hoxtnyics a> imunnwlure antvbiianshi es THECI » 











one 


4th card: 


VKTS 
FORMAT (F 10. 2) 
VKTS is velocity of boat in knots. One card is required for 


each speed. Number of cards equals NUMBER on card 3, 


afl 


CARD 2 


CARD 1 


BINARY DECK 






+ BINARY 


Fi gure 23, Program Assembly. 
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STABILITY CHARACTERISTICS FOR PLANING BOAT SERIAL 6809 


EQUILIBRIUM PLANING CONDITIONS 


VELOCITY (FPS) = 25, 33 

TRIM (DEG. ) = 4.275 

ASPECT RATIO = 4.03 

WETTED KEEL (FT) = 5.02 WETTED CHINE (FT) = 4.58 
WETTED AREA (FT**2) = 5,73 DRAG (LBS) = 21.40 

DRAFT (FT) = .37 

MEAN WETTED LENGTH (FT) = 4. 80 


ESTIMATED EHP = . 99 


STABILITY INDICATOR = .18479E-01 


Figure 25. Sample of Computer Output when both LIST 1 and 
LIST 2 equal 2. 


000 J/ISEE TASH DIN SP 2< HOY P°PYTERI@T AHS YTLAIGATC | 










PACITY I OMANI MULLS 


ee 88 = (e@D YTMOUaV 
£0.b > OIPAR TSUG 
| $0.8 = (T4) ey GaPTa 
wARG tr .8 = (S8F%) AgaA GAT 
TS. = (PE) 7 
98.) = (TS) ROMA CATTOwW BA! 


. 








ye. & "ae" 


)= 
Py =: 


eas = (aul) 


fe ( 4 









= = 


* 








206 Ons 


COMPUTER PROGRAM TwO 


STAB Y OF PLANING sGRaiae, 

DIMENSION T1(15)»s VALUBMB yes IDENT (5) 

COMMON AAs Alls A225 ASP» BBs Bll» B22s BEAMs BETAs CC» 
Melis @22s CL Bs CLO» CGs 9CVemDDs Clls D22s)DifBis DIEGis 
ZIeRAGs EEs Ells E22» EPSILs Fs Gll» G22s LISTls LISTI2s NR 
3RHO» Se» TAUs TRIM» Tl» VALUEs Vs VCGs» VKTS» VMs Ws WETL» 
4WCHINEs WKEEL,» YI 


READ lsLIST1ls LIST2»NBOATS 


BewEISTI= 1» PRINT OUT ALL COEFFICIENTS AND DERIVITIVES 
ASSOCIATED WITH STABILITY EQUATIONS 

IF LIST 1 = 29 PRINT OUT ONLY STABILITY INDICATOR 

IF LIST ¢ = Ll» PRINT OUT ALL ITERATIONS INVOLVED IN 
SOLVING FOR EQUILIBRIUM PLANING CONDITIONS 


IF LIST 2 = 29 PRINT OUT ONLY FINAL PLANING CONDITIONS 


FORMAT (313) 

DO 7 NN=l»NBOATS 

READ2 sBETAI »sEPSILI sF 9 VCG BEAM sCGoRHOsYI oW 
FORMAT(4F50295F 1002) 

BETA=BETAI/57e2956 

EPSIL = EPSILI/57e2956 

READ 3»sNUMBER»s (IDENT(1I)sT=195) 

FORMAT (1395A4) 


PRINT1IO1ls(IDENT(I)»sT=195) 





a7 jee 


DO 7 M = 1»sNUMBER 
READ 45 VKTS 
Vo = VKTS* 14689 

4 FORMAT(F10e2) 
CALL ANGLES 
PALL COEFF] 
EHP=V*¥DRAG/550 « 
PRN le hG EHP 
RD=BB*¥CCHEDD—-AA#D0D*#*2-BB**2*EE 
ieeGtel- Si llores oO 9 © 

> PRINT102sAA sBBeCCsDDsEE 
PRINT1039A11 
PRINT1049B11 
PRINT1059C11 
PRINT1069D11 
PRINT107sE11 
PRINT1089611 
PRINT1099A22 
PRINT1109B22 
PRINT111sC22 
PRINT112sD22 
PRINT TYom2 22 
PRINT114sG22 

SePRINTII5 »RD 

7 CONTINUE 

101 FORMAT(43H1ISTABILITY CHARACTERISTICS FOR PLANING BOAT» 


19H SERIAL »5A4///) 


HOZePOKMAT(4S6H THE COEFFICIENTS OF THE FOURTH-ORDER EQUATIONS 





ete 


153H*OF MOTTON 


2933H 060) ARE AS FOLLOWS (AA THRU 

mos FORMAT(G6H FORCE DERIVIMMEVE FOR A 
LISHICAL ACCELERATION =9F10.3) 

104 FORMAT(46H FORCE DERIVITIVE FOR 

LIOHICAL VELOCITY =,F 10.3) 


105 FORMAT(46H FORCE DERIVITIVE FOR 


LIOHICAL POSTTION »9F 10.3) 


106 FORMAT(46H FORCE DERIVITIVE FOR 
L1OHLAR ACCELERATION =sF10¢3) 
107 FORMAT(46H FORCE DERIVITIVE FOR 
L1OHLAR VELOCITY =9F 106.3) 
108 FORMAT(46H FORCE DERIVITIVE FOR 
L19HLAR POSITION =9F10.3) 
109 FORMAT(47H MOMENT DERIVITIVE FOR A 
L1OHLAR ACCELERATION =9F10.}) 
110 FORMAT(47H MOMENT DERIVITIVE FOR 
119HLAR VELOCITY =3F10.3) 
111 FORMAT(47H MOMENT DERIVITIVE FOR 
L1OHLAR POSITION =9F10.3) 
112 FORMAT(47H MOMENT DERIVITIVE FOR 
LIOGHICAL ACCELERATION =9F10¢3) 
113 FORMAT(47H MOMENT DERIVITIVE FOR 
LLOHICAL VFLOCITY | =9F 10.3) 
114 FORMAT(47H MOMENT DERIVITIVE FOR 
LIOHICAL POSITION =9F 10.3) 
115 FORMAT(22H STABILITY INDICATOR 


116 FORMAT(1€4 ESTIMATED EHP =sF20e2// 


EEN SE 2564/7) 


UNTT 


UNTT 


UNIT 


UNTT 


UNIT 


UNIT 


ENGREMENT 


INCREMENT 


INCREMENT 


INCREMENT 


INGREMENT 


INCREMENT 


(AA*S**4 + BER*S*#*S + GCENSEAZ + DDORS + 


OF 


OF 


OF 


OF 


OF 


OF 


UNIT INCREMENT OF 


UNIT 


UNIT 


UNIT 


UNIT 


UNIT 


/ ) 


INCREMENT 


INCREMENT 


PNG REMEND 


INCREMENT 


INCREMENT 


=sE15,.5////) 


OF 


OF 


Or 


OF 


OF 


ae 


VERT » 
VERT » 
VERT» 
ANGU » 
ANGU » 


ANGU » 


ANGU >» 


ANGU » 


ANGU » 


VERT» 


VERT » 


VERT » 





CALE Exit 


END 











ital 


SN oes 


SUBROUTIN= ANGLES 


DIMENSTON 71415) »VALUB Gi) 

COMMON AAs Alls A22sx ASP» BBs f}1ls B22x BEAMs BET As CC» 
1Glis ‘C225 CLBs CLO>s CGgmeG@y DD, Dilys sD22,) Dib bic 
ZDORAGS EEGs Ells F22s EPSTEsg=s Glls G22s 5 UISils LISiZzs ans 
3RHO» S» TAUs TRIMs Tls VALUEs Vs VCGs VKTS»s VMs Ws WETL» 
4SWCHINE» WKEEL»s Yi 

CV=V/SQRTF(3262*BEAM) 

CLB=2e*W/ ( RHO*( (V¥BEAM) #¥*2 ) ) 

BETAI=BETA*57.22956 

ICLO=1 

eEOl={1¢ELloO 

CLO=CLO1*0.201 

IF (CLO-O0 -O065*BETAI#¥CLO#*.6-CLB) 29394 

heme =1C LOD 

eo TO 5 

CLO=CLO-.001 

BEC CEO—05006 S8BETAI*¥CLO®*.6-CLB) 35356 

GO TO 4 

CLO=CLO 

N=1 

NR = N 

TLIN)=le 

PEO hSWZel)lisisll 

PRINT 1000sCLB»CLOsT1(N) »N 

CALL FACTOR 


IF (ABSF(VALUE(N) )-20001)1091097 





Tae 


TLIN)=200 
(PULNST 21 | Aslan 
13 PRINT 1000%CLBsCLOsT1IN) oN 
MP CALL FACTOR 
IF (ABSF (\ ALUE(N) )—00001)1091098 
8 CONTINUE 


DO 9 N = 3598 


NR = N 
SLOPE = (VALUE(N-1)}) -— VALUE(N-2))/0T1°N-1) =-T1I(N=2)) 
ic N) = TitUN=))— VALUE(NS) 7 SeOre 


[ecclSl2=1) 15914515 
14 PRINT 1000sCLBsCLOsT1(N)3N 
15 CALL FACTOR 
IF (ABSF(VALUE(N) )—20001)1051099 
9 CONTINUE 
10 TERM=e5*SINF( BETA) *COSF(TAU)/3¢1416/S INF (TAU) /COSF( BETA) 
WKEEL=BEAM* (ASP+TERM) 
WCHINE=BEAM* (ASP=TERM) 
DRAFT=WKEEL¥SINF( TAU) e 
WETL=ASP*BEAM 
TRIM=TAU*5722956 
S=WETL¥BEAM/COSF (BETA) 
PRINT 1002 
PRINT 1003 9VsTRIMsASP sWKEEL sWCHINEsSyDRAGsDRAFT sWETL 
1000 FORMAT( 26H CALLING FACTOR WITH CLB =9F9e493Xs7H »sCLO =s 


1F9e¢493Xs13H »sAND TRIM =9F2004515X94H N =913) 





1002 FORMAT( 31H 
1003 FORMAT (17H 
11511 ASPECT 


2200 WETTED 


ec 


EQUILIBRIUM PLAN! NG CONDITIONS//) 
VELOCITY (FPSie=5F 7.2/1 4R RIM E Ge) sour 
RATIO =sF5e27190H WETTED KEEL (FT) seg Fos2s207 


CHINE (FT) =5F66.2/22H WETTED AREA (Flex jenn 


SF ie2sl2Xe1l3H DRAG (LBS) =sF7e2/13H DRAFT (FT) “2 sFS5a27 


426H MEAN WETTED LENGTH (FT) =s9F6e0e2///) 


RETURN 


END 





= 


SUBROUTINE FACTOR 


DIMENSION T1i(15)5 VALUES 
COMMON AAs Alls A22s ASP» BBs Blls B22s BEAMs BETAs CC» 
1Ells G22s CLBs CLOs CGs EVs DDs Dlls D225 Dito) ss CIFGiE 
ZORAGs EEs Ells E22s EPSIEsats Gll»s G22s Blois CIS! Zeus 
3RHOs Ss TAUs TRIMs Tls VALUTs Vs VCGs VKTSs VM» Ws WETL» 
GWCHINEs WKEEL»s YI 
N = NR 
IF (T1(N))29293 
2 ITF(LIST2-1)2294022 
4 PRINT 1009NoT1(N) 
@eetliN) = «5 
3 ASP = 80e*(CLO/T1(N)*¥*¥1.1) 
TAU=T1(N)/57e2956 
IF ( ( ¢O12*SQRTF (ASP) +e0055#ASPHH2.5/CVHEH2Z)*XTLIN) ¥¥#1 6-1-CLO) 
1118512 
11 ASP=ASP+,01 
IF ( (20123 SQRTF (ASP) +e0055*# ASP##2 -5S/CV¥¥2)*¥#TI(N) ¥¥16¢1-CLO) 
111898 
iemasoP = ASP - .01 
IF ( (eO012*SQRTF (ASP) + e005 5*# ASP #2 5 /CVHEK2)*¥TI(N) ¥*#1¢1-CLO) 
1858912 
B C=CG—( eo 75-le/(5a21* ( CV/ASP ) #¥ 242639) ) ¥ASP*#BEAM 
VM = V¥SQRTF (1lemeO1lL2Z*TAU*#1] CL /SQRTF (ASP) *(85e-50e* BETA) / 
1COSF (BETA) #*2/COSF( TAU) ) 
RE=131770e* ASP*BEAM*V 


REE=e43429448*LOGF (RE) 





50 


ole, 


60 


65 


70 


1, 


80 


85 


90 


95) 


ae fee 


CF=.008179 

DEL=.001 

CF=CF=}{DEL 

FRE=e¢242/SQRTF (CF )-24329448*LOGF (CF) 

CGR Bik E555 951s 50 

ITF (DEL—e001)65 360960 

CF=CF+DEL 

DEL=e0001 

GO TO 50 

fiterwDelb=s00C1) 75570970 

CF=CF+DEL 

DEL=.0000) 

GO TO 50 

IF (DEL—e00001)85s80s80 

CFaCF+DEL 

DEL=.e000001 

Ge TO 50 

IF (DEL—2e000001 )95390990 

CF=CF+DEL 

DEL=e0000001 

GO TO 50 

CFT=CF+e0004 

A=VCG~-e25*BEAM*SINF(BETA)/COSF (BETA) 

DRAG =RHO*(VM#BEAM) #*2#ASP#CFT#.5/COSFI(BETA)/COSFI(TAU) + 
LW*SINF(T/ J) /COSF (TAU) 

VALUE (N) =W*¥(C¥(COSF(CEPSTL)=SINF (TAU) #SINF(TAUtEPSIL) )-F* 
PSINF (TAUJMCOSF (TAU) )4DRAGEICH(SINE (TAUREGOSF(EPSILE}= 


2SINF(TAU+EPSIL) )+(A¥COSF (CEPSIL)—-F)*COSFI( TAU) ) 





Eo foame 


LRU sh S a 2 li) Ges 916 
5 PRINT 101 sVALUE(N) sCFT 
100 FORMAT(47H NEGATIVE ANGLE ENCOUNTERED ON ITERATION NUMBER, 
113%940H THE VALUE OF TATSBANGLE (INSDEGRERe) = 1SsFl0se) 
101 FORMAT(32H RETURNING TO ANGLES WITH VALUE#sE12¢5>» 
PrOnPAND® CET =sET2.5/777) 
6 RETURN 


END 





Gy) (3 


a@ 


SET Cie 


SUBROUTINE wGG@ Gira 


DIMENSION T1(15)» VALUE(15) 

COMMON AAs All» A229 ASP» BB» Bll» B22» BEAMs BETAs CC» 
MNGi) C229 LBs CLOs COMMGUREDDs DilelD2 a eDiime ie Dneoie 
PRAGHMEEs Ell, E22s EPSTMENE Glls Goon TncimenGT NMA 
BRHGbmeSs TAU, TRIMs Tis) VAlUem Vs VGGy ViKTSiuyemnimmncn nee 


GWCHINEs WKEELs YI 


EQUATIONS OF MOTION OF PLANING HULLS 
AA»sBBsETC=ROUTH CRITERION FACTORS 
Al1sB11lsETC=NONDIMENSIONAL FORCE AND MOMENT COEFFICIENTS 


AlsBlsETC=FORCE AND MOMENT COEFFICIENTS 


EPS = EPSIL 

A = ASP 

VERM=0125*RHO*3 o1416%(WETL*¥#*2) #BEAM 

VERY 1 =00625%* 60625*RHO#3 0141 6% ( WETL*¥*4) *BEAM 


CPLEWETL*#(00¢75-1e0/ (5021 *( (CVA) #¥2) +2039) ) 


Pe PeERENTIATION OF LIFT COEFFIGTENTIMXSAREA WITH RESPECT 


TO TAU 


TAUL=TAU-0.2001 


TAU2=TAU+0.001 


WC WCHINE 
WK = WKEEL 


A = 1e/ASP 








- 80 - 


CLVOLIl=Hle/(2e*WETL¥(CVHE2) ) *( (CC WCHEH2)HSINE( 2e*TAUL) 
L/BEAM+1¢/30*(2e¢*WC+WK) *SINF( BETA) /COSF (BETA) )) 
BS GEB1l=065¥ELV0L1 
CLSA1=(1e35708*A*TAUL* (COSF ( TAU] )*#2)#(1.0—-SINF( BETA) ) 
1/(1¢0+A)+4e0*(SINF( TAU) ¥*2)*(COSF(TAUL)**¥3)*COSF(BETA)/ 
Z25<e0tCLB1)*3S 
CLVOL2Z£le/(2e*WETL¥(CV##2) )#( CC CWCHH2) SINE (20* TAL 2) 
L/BEAM+le/3e*(2e*WCH+WK) *¥SINF (BETA) /COSFI(BETA))) 
peCLB2=0e5*CLVOL2 
CLSAZ=(1-¢5708*¥ A*TAUZ*(COSF (TAU2) #¥*2)*(1e¢0-SINFI(BETA) ) 
1/(1e0+A)+4e0*(SINF (TAUZ) #2) *( COSF(TAU2)*¥*3)*COSF( BETA) / 
woe OFCl a2) *S 
CL=1.5708* TAUXA*#COSF (TAU) ##2*#(1e-SINFIBETA) )/7(1letA)+4e% 
LSINFC TAU) **®2*COSF( TAU) ##3*COSF(BETA)/3e+(CLVOLI+CLVOL2 )/4e 
TAO=TAU1 
TABHETAU2 
CO=1¢5708* TAUXA#COSF (TAO) ##2*(Le-SINF( BETA) )/7(1LletA)+4e* 
LSINF (TAO) **®2*COSF (TAO) ##3*COSF (BETA) /36¢+(CLVOL1)/26 
CB=1¢5708* TAUXA*COSF (TAB) #¥#2*(Le-SINFIBETA) )/(1letA)+4e* 
LSINF (TAB) ¥*2*COSF (TAB) **#3*¥COSF (SETA)/3e+(CLVOL2Z)/26 
DADTAU=1000 e#W/ (RHO*#¥V##2 )#(16/CB-1¢/CO) 


DIFB1=(CLSA2Z2=CLSA1)/e002+CL¥DADT AU 


DIGPPERENMeyATION OF LIFT CORPSE Awith RESPECT [O° Z 


DELZ=0.001 
DELA=-DELZ*BEAM/SINF (TAU) /WETL¥*2 


DELS=BEAM*DELZ/(SINF( TAU) *COSF (BETA) ) 





ee be 


DEEE=DElz SINE (Au 
CLVOL3=( ((WC-DELL) #*2) *SINF (TAU) /BCAM4+(26*WC+WK=3 6 *DELL ) 
L*SINF (BETA) /(COSF(BETA)*3.))/(2e*(WETL=DELL) *(CV#¥#2) ) 
We CHB = 0S %GiEVOL3 
CLSA3=(1+65708* (A=DELA) *TAU#(COSF (TAU) **#2)*(1e-SINFI(BETA) ) 
L/(LetA-DELA) +4 e#(SINF (TAU) #¥#2) *(COSF( TAU) **3) *#COSFIBETA)/3.¢ 
2+CLB3)*(S-DELS) 
CLVOL4=(( (WC+DELL) **2) *SINF (TAU) /BEAM+(20¢*WC+WK+3 6 *DELL) 
L¥SINF (BETA) /(COSF (BETA) #36))/(2e*(WETL+DELL) *(CV¥*2) ) 
9 CLB4=0254 ZLVOL4 
CLSA4=(125708* (A+DELA) *TAU®(COSF (TAU) **2)*(1e-SINF( BETA) ) 
L/CLetA+DELA) +4o%(SINF (TAU) #¥2) *(COSF (TAU) **3)*COSF (BETA) /36 
2+CLB4)*(S+DELS) 


DIFCl=(CLSA4=-CLSA3)/0.002 
BTrrFERENTIATION OF MOMENT WITH RESPECT TO Z 


CPL1=075-(WETL—DELL ) /(5621*(CV*BEAM) **2/ (WETL-DELL) *#2+2 039) 
CPL2=.75—-(WETL+DELL ) / (521% (CV¥BEAM) ¥#2/ (WETL+DELL ) ¥¥2+2 039) 
DCPLDZ=(CPL2=CPL1)/.002 

C1=065*RHO*(V¥*2)*OIFC] 
DMDZ=(CG=CPL) #C1/COSF (TAU) -W#¥DCPLDZ/COSF (TAU) 

Al=W/3262+VERM 

B1=0+5*RHOXV*DIFB1 

C1=0.5*RHO* (V¥*2)*DIFC] 

D1 =VERM¥VCG* (COSF ( TAU) ¥(CG—.5*WETL) /VCG-SINF( TAU) ) 
E1©B1*(COSF( TAU)* (CG-WETL/20)-VCG*SINF (TAU) } 


Gl=V*B1 





= uy es 


A2=VI+VERYI 

B2=E1*( (CG—<CPL) *COSF ( TAU)=—VCG*¥SINF( TAU) ) 

C2=G1¥*¥( (CG=-$CPL )*COSF(TAUI-VCG¥SINF (TAU) )+WH#(CPL¥(SINF (TAU) =VCG* 
LCOSF (TAU)/SSINFE (TAU) /WK#(COS*(TAU)/SSINFE (TAU) =1Le/SINF( TAU) ) )=-CG* 
Z2SINF(TAU)-VCG*¥COSF (TAU) ) 

D22VERM* (COSF( TAU) *(CG=WETL/2¢ )-VCG*SINF (TAU) ) 

e2=B ret (CG—CPL *COSF (TAU eV GGrS TINE TAU) 

G2 = DMDZ 

El=E1+V*VERM 

G1=G1+V*81 

B2=B2+V*D2 


C2=C2+V¥E? 


U IS AN ARBITRARY NONDIMENSIONALIZING VELOCITY 


YU = 10e 
QUANTITY O25 RHO CANCELLED OUT OF ALL FOLLOWING 


Al1l=A1/(BEAM#*3 ) 

Pl = 617 (Ur BEAMAE2 )) 
CLI=HEC1L/(BEAM* (U##2) ) 
D11=D1/(BEAM##4) 
ELL=E1/(U*(BEAM**3 ) ) 
G11=G1/ ( (BEAM*U) ##2) 
A22=A2/(BEAM**5 ) 
B22=B2/(U* (BEAM##G) ) 
C22=C2/( (BEAM*##3 ) % (U*#2 ) ) 


D22=D2/(BEAM##4) 





Seto 


E22=£2/(U*( BEAM*#3 ) ) 

G22=G2/( (BEAM*U) **2) 

AA=le 
BB=(A22*B11+A11#*B822-D22*E11-E22*D11)/(AL1*A22-D11*D22) 
CC=(A22*C114+B22*B11+A11¥*C22-D22*G11-E22*E11-G22*D11)/(Al1 
1*A22—-0D11*D22) 
DD=(822*C11+B11*C22-£22*G11-G22*E11)/(A11*A22—-D11*D22) 
EE=(C22*C11-G22*G11)/(A11*A22-D11*D22) 

RETURN 


END 
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